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Abstract. Neutron star binaries, such as the one observed in the famous 
binary pulsar PSR 1916+13, end their life in a catastrophic merge event 
(denoted here NS^M). The merger releases ~ 5 • 10^"^ergs, mostly as neutri- 
nos and gravitational radiation. A small fraction of this energy suffices to 
power 7-ray bursts (GRBs) at cosmological distances. Cosmological GRBs 
must pass, however, an optically thick fireball phase and the observed 7-rays 
emerge only at the end of this phase. Hence, it is difficult to determine the 
nature of the source from present observations (the agreement between the 
rates of GRBs and NS^Ms being only an indirect evidence for this model). 
In the future a coinciding detection of a GRB and a gravitational radiation 
signal could confirm this model. 



1. Introduction 

Binary neutron star merger (NS^M) is the last event in the complex evolu- 
tion of some massive binary systems. A massive binary becomes a massive 
X-ray binary after one of the stars undergoes a core collapse and a su- 
pernova explosion. The second supernova disrupts most systems. A small 
fraction survives and forms a neutron star binary. There is a chance to 
detect a neutron star binary if one of the neutron stars is a pulsar. So far 
three such binaries have been detected. General Relativity predicts that 
the binary will emit gravitational radiation and its orbit will shrink. Taylor 
and Weisberg [1] measured the orbit of PSR 1916-|-13 and confirmed this 
prediction with an amazing precision. The gravitational radiation emission 
from the known binary pulsars is too low for a direct detection. The emis- 
sion increases as the distance between the two neutron star decreases and 
gravitational radiation detectors such as LIGO [2] and VIRGO [3], that are 
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being built now, will be able to detect gravitational signals from the last 
three minutes [4] of binary neutron stars from cosmological distances. 

The outcome of the merger is most likely a rotating black hole of ~ 
2.2 — 2.8m0 (a massive rapidly rotating neutron star, very near the upper 
mass limit for a rotating neutron star, cannot be ruled out yet [5]). In 
either case the released binding energy is ~ 5 • 10^"^ ergs. Shortly after 
the discovery of the first binary pulsar, Clark and Eardley [6] estimated 
that most of this energy is emitted as neutrinos, and a smaller fraction as 
gravitational radiation. The neutrino signal resembles a supernova neutrino 
burst. At present there is no way to detect such bursts from cosmological 
distances. Even if such a detector would have existed it would have been 
impossible to distinguish the rare NS^Ms neutrinos bursts from the much 
more frequent SN neutrino bursts. 

NS^Ms conspire to release their energy in the two channels that are hard- 
est to detect: gravitational radiation and neutrinos. However, if a small frac- 
tion of the total energy is channeled into electromagnetic energy it would 
be detectable. Several years ago Eichler, Livio, Piran and Schramm [7] sug- 
gested that vv annihilation {vv e'^e~) can convert ~ 10""^ of the total 
neutrino energy to electromagnetic energy which in turn would produce 
a GRB (a possibility mentioned without a detailed discussion earlier by 
Goodman [8] Paczyhski [9] and Goodman Dar and Nussinov [10]). Accord- 
ing to this model GRBs signal the final stage in the complicated evolution 
of what was initially a massive binary system: the merger of a neutron star 
binary into a black hole (or a rapidly rotating maximal neutron star). 

This idea was accepted with skepticism. At the time it was generally 
believed that GRBs originate from neutron star in the Galaxy (see e.g. a re- 
view by Higedron and Lingenfelter [11]). However, the BATSE detector on 
the COMPTON-GRO observatory have revolutionized our understanding 
of GRBs [12, 13]. The BATSE observations have shown that the angular 
distribution of GRBs is isotropic and the count distribution is incompatible 
with a homogeneous one. The isotropy and inhomogeneity rule out a Galac- 
tic disk population, which would be either anisotropic (if the sources are at 
typical distances of more than several hundred pc) or homogeneous (if the 
sources are at typical distances of less than several hundred pc). This leaves 
us with the possibility of a Galactic halo population (with a very large core 
radius) or a Cosmological population. While the controversy between a 
Galactic (mostly Galactic halo) and Cosmological populations is still going 
on, there is a growing consensus (see e.g. Fishman [13]) that GRBs are 
cosmological. With this emerging consensus the NS^M model has made the 
full transition from one of the least likely amongst more than a hundred 
GRB models to the most conservative one. It is the only model based on 
an independently observed phenomenon which, as we discuss shortly, takes 
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place at a comparable rate and can easily account for the energy required. 

In this lecture I discuss three issue. I examine first the GRB distribution 
and I show that the count distribution is compatible with a cosmological 
distribution (an arbitrary inhomogeneous distribution will not necessarily 
be compatible with a cosmological one). The analysis of the count distri- 
bution enables us to estimate the rate of GRBs, which I compare with 
the estimated rate of NS^Ms. The agreement between the two rates clearly 
supports the NS^M model. 

I turn, next, to the issue of fireballs: The sudden release of copious 7- 
ray photons into a compact region, as required for a cosmological GRB, 
creates an opaque photon-lepton fireball due to the prolific production of 
electron-positron pairs. Once this was an argument against any cosmolog- 
ical source[14, 15]. However, Goodman [8] and Paczyhski [9] have shown 
that the fireball expands relativistically and it releases its energy at a lat- 
ter stage when it is optically thin. This optically thick phase between the 
energy source and the final emission stage makes the task of deciphering 
the GRB enigma much more difficult that what was otherwise expected. 

Finally, I turn back to NS^Ms. I discuss some recent numerical simu- 
lation of the hydrodynamics of the merger and I examine how does this 
model satisfy the constraints introduces earlier. I conclude with some open 
questions and predictions. 

2. The Distribution of GRBs 

2.1. ANGULAR DISTRIBUTION 

The angular distribution of the bursts is isotropic to within the statistical 
errors of the sample [13]. This is clearly in agreement with a cosmological 
model. The isotropy sets some sever constraints on Galactic halo models 
[16] requiring a large homogeneous core (significantly larger than the ho- 
mogeneous core expected in the dark matter distribution) and pushing the 
GRBs to almost inter-galactic distances. The lack of observations of GRBs 
from M31 sets an upper limit to the possible distance to GRBs and the 
combined observations begin to rule out Galactic halo models (see however 
[17] for a recent discussion of Galactic models). 

2.2. COUNT DISTRIBUTION 

A homogeneous population of sources in an Eucleadian space-time will have 
a count distribution N(C) oc C~"^/^. Several relativistic effects infiuence a 
cosmological distribution of sources. First, space is not Eucleadian and the 
relation F oc which is the basis for N(C) oc C~"^/^ is not valid. Second, 
the redshift factor, 1 + z, changes the observed spectrum and introduces a 
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K correction. Third, the number of counts is diluted by a (1 + z) factor if 
the detector operates within a fixed interval AT that is shorter than the 
total duration of the burst. These effects combine to yield: 

for a detector with a fixed energy range, A_E that operates for a fixed time 
interval, Ai, and sources with a count spectrum: N{E) oc E~" . L depends 
on the luminosity of the source, i, in the relevant energy range, A_E, on 
the average energy E and on the observation time, Ai: L = L(AE)At/ E. 
di(z) is the luminosity distance [18]. The number, N(> C), of events with 
an observed count rate larger than C is: 

foo fz(C,L) J2 J„ / ^\ 

N(>C) = 4.l n(L)4Ll ^n(,)^,„ (2) 

where rp(z) is the proper distance to a redshift z and z(C,E) is obtained 
by inverting Eq. 1. n(E) is the luminosity function and n(z) is the intrinsic 
rate: the number of events per unit proper volume and unit proper time at 
redshift, z. N(C) depends on the cosmological parameters: i7 and A, and 
on the source parameters: a, n(z) and n{L). 

Cohen and Piran [19] calculated the likelihood function that the BATSE 
2B (the 2B catalogue) data results from a cosmological distribution of the 
form given by Eq. 2 for a variety of cosmological models and source param- 
eters. Following Kouveliotou et. al. [20] and Mao, Narayan and Piran [21], 
Cohen and Piran divided the GRB population to two sub -populations, 
short (^^90 < 2 s) and long (^igo > 2 s) bursts (BATSE is more sensitive to 
the long bursts [21] hence it is meaningless to perform the analysis on the 
whole population as one group). 

The maximal red shift up to which BATSE detects long bursts (for stan- 
dard candles with no source evolution and spectral index a = 1.5 [22]), is 
Zmax{long) = 2.11q'7. With estimated BATSE detection efficiency of ~ 0.3 
this corresponds to 2.3^°']' • 10~^ events per galaxy per year (for a galaxy 
density of 10~^h'^ Mpc""^ [23]) the rate per galaxy is independent of Hq 
and it is only weakly dependent on i7. For J) = 1 and A = the typical 
energy of a burst whose observed fiuence is Fr (in units of 10~ ''ergs/cm^) 
IS These numbers vary slightly if the bursts have a wide 

luminosity function. The distance to the sources decreases and correspond- 
ingly the rate increases and the energy decreases if the spectral index is 2 
and not 1.5. The rate increases and the luminosity decreases if there is a 
positive evolution of the rate of bursts with cosmological time. 

Short bursts are detected only up to a much nearer distances: Zmaxishort) 
= 0.4"'"^'^, again assuming standard candles and no source evolution. There 
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is no significant lower limit on z^ax for this sub-class. The estimate of 
Zmaxishort) Corresponds to a comparable rate of 6.3~^'^ • 10~^ events per 
year per galaxy and a typical energy of 3"*""^^ • 10'*^i^_7 ergs (note that there 
is no lower limit on the energy or upper limit on the rate since there is no 
lower limit on Zmaxishort). 

Luminosity functions with an effective width of up to factor of ten fit 
the data. This width is comparable to, and even wider than, the width of 
some observed luminosity distributions such as the luminosity functions of 
different types of supernovae. 

Norris et. al. [24] found that the dimmest bursts are longer by a fac- 
tor of ~ 2.3 compared to the bright ones With our canonical value of 
z-maxi^ong) = 2.1 the bright bursts originate at Zfj^ighti^ong) ~ 0.2. The 
corresponding expected ratio due to cosmological time dilation, 2.6, agrees 
with this measurement. If the interpretation is correct than this result con- 
firms the cosmological model. It implies that intrinsic source evolution is 
insignificant and it rules out the low z^^x obtained form the combined PVO 
and BATSE data [25]. 

3. Fireballs 

A cosmological GRB releases lO^^ergs (if the energy emission is isotropic) 
in a very small volume. The rapid rise time observed in some bursts implies 
that the sources are compact with sizes, Ri, as small as 100km. This results 
in what we call a "fireball" : an optically thick radiation - electron - positron 
plasma whose initial energy is larger than its rest mass. The initial optical 
depth in cosmological GRBs for 77 e~^e~ is: 

T^^ = fgEcjTlR^m.c^ ^ lO^'^fgE.^siR-h (3) 

where Ei^^i is the initial energy of the burst in units of lO^^ergs, Rij is the 
initial radius in units of 10''cm and is the fraction of primary photons 
with energy larger than 2meC^. Since t^^ ^ 1 the system reaches rapidly 

thermal equilibrium with a temperature: T = QAE^^'^R- ^''^MeV. At this 
temperature there is a copious number of e"*" — e~ pairs that contribute 
to the opacity via Compton scattering. (It is interesting that Eq. 3 yields 
T^-y ^ 1 even for Galactic halo objects [26]). 

The huge initial optical depth prevent us from observing directly the 
radiation released by the source. The observed 7-rays emerge only after 
the fireball has expanded significantly and became optically thin. The fire- 
ball phase determines, therefore, the observational features of GRBs and it 
screens the specific nature of the energy source. A comparable situation is 
familiar in stars where the energy generated in the core leaks out through 
an optically thick envelope and the observed spectra is independent of the 
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details of the energy generation mechanism at the core. Another analogous 
situation occurs in SNRs where the observed radiation depends just on the 
total energy of the ejected material and is independent of any other specific 
features of the source. I review here some essential physics of fireballs (see 
[27] for details). 

3.1. FIREBALL EVOLUTION 

Consider, first, a pure radiation fireball. Initially, the local temperature 
T ^ nieC^ and the opacity due to e+e" pairs, ^ 1 [8]. The fireball 
expands and cools until at T,p ~ 20 KeV, ~ 1 and the photons escape 
freely as the fireball becomes transparent. 

Astrophysical fireballs include baryonic matter in addition to radiation 
and e+e" pairs. The baryons affect the fireball in two ways: The electrons 
associated with the baryons increase the opacity, r = + r^, delaying the 
escape of radiation. The baryons are also dragged by the accelerated leptons 
and this requires a conversion of the radiation energy into a kinetic energy. 
Thus, two important transitions take place as a loaded fireball evolves: The 
transition from optically thin to optically thick fireball, which takes place 
at: 

i?. = (^)^/^ = 6.10-cmO"^ (4) 

and the transition from a radiation dominated phase to a matter dominated 
phase which takes place at: 

Rr, = 2R,7] = 2 X lO^^cm R^jT]^. (5) 

T] = Ei/Mc^ , the ratio of the initial radiation energy E to the rest energy 
Mc^ controls these transitions. 

The overall outcome of the fireball depends critically on whether R,^ > 
Rt or vice versa. If R,^ > i?^, most of the energy comes out as high 
energy radiation, while if R^- > i?^, the fireball results in a relativis- 
tic expanding shell of baryons. Energy conservation dictates that in this 
case M^jp ~ E and ^jp ~ r]. We can classify four situations [28, 29]: (i) 

1^ > r]pa^r = (Sa^E.aT.^ / AiTmlc'^ R. f ^ IO^^eHIr'^^^: The effect of the 
baryons is negligible. The pair opacity Tp drops to 1 and ^ 1 at Tp while 
the fireball is still radiation dominated and the radiation escapes carrying 
all the energy, (ii) 1]pa^r > V > Vb = {SaTE^STrnipC^ Rf^/^ ^ W^eIHr-^^^: 
The opacity becomes dominated by r^. The comoving temperature de- 
creases far below Tp before r reaches unity. The fireball is, however, still 
radiation dominated when = 1 and the escaping radiation carries most of 
the energy, (iii) r]b > r] > 1: The fireball becomes matter dominated before 
it becomes optically thin. The total energy is the bulk kinetic energy of 
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extreme relativistic baryons. The final Lorentz factor is ^jp ~ r]. (iv) 77 < 1: 
This is the Newtonian regime. The rest energy exceeds the radiation en- 
ergy and the expansion is not relativistic. This is the situation, for example 
in supernova explosions in which the energy is deposited into a massive 
envelope. 

A quick glance at the corresponding mass limits (the transition from 
case (ii) to case (iii) for Ei ~ lO^^ergs is for M < 1Q~^Mq) reveals that 
case (iii) is the most likely one and even this requires a rather "clean" 
fireball (The transition to the Newtonian regime is at M = 3 • IQ-^Mq). 
Initially such a fireball is radiation dominated and it accelerates with 7 oc r. 
The fireball is roughly homogeneous in its local rest frame but due to the 
Lorentz contraction its width in the observer frame is Ar ~ Ri, the initial 
size of the fireball. From R,^ onwards the baryons coasts asymptotically 
with 7j? ~ T]. 

3.2. ENERGY CONVERSION MECHANISMS 

The kinetic energy of the bulk motion of the relativistic baryons can be 
recovered if shock waves form. The shocks could be either internal [30, 31, 
32] or due to interaction with the ISM [33, 34, 35]. Variations of r] (and 
hence in 7^?) as a function of radius become important at i?^,: 

R^ ^ r]^R, ^ lO^^cm 7. (6) 

If 7j? decreases outward than inner shells take over outer shells and internal 
shocks form [30]. Quite generally this is preceded by an unstable phase 
[32]. The shocks [30, 31] or the instability that preceded them [32] could 
convert a significant fraction of the kinetic energy of the baryons back to 
thermal energy. Since Ryj > R^- (for most reasonable parameters) this takes 
place in an optically thin region and the emitted photons can escape freely, 
producing the observed GRB. 

If 7j? increases monotonically outward than the observed pulse width 
increases linearly with the radius: Ar ~ Ri/^jp and there are no internal 
shock. In this case the baryons can still interact with the ISM [33, 34, 35]. 
Similar situation occurs in SNRs where the interaction of the ejecta with the 
ISM produces a shock in which the kinetic energy of the ejecta is converted 
into radio emission. The mean free path of a relativistic baryon in the 
ISM is ~ lO^^cm, hence the interaction between the baryons and the ISM 
cannot be coUisional. However, from the existence of SNRs we can infer 
that a coUisionless shock can form (possibly via magnetic interaction). The 
interaction becomes significant at R^ where the fireball sweeps an external 
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mass of Mq/jp = (Ei/r]c'^)/jp and looses half of its initial momentum: 



Mo 



1.3 X lO^^cm E-'^-{r]^^'-^n 



(47r/3)ra7i?. 



l/3„-2/3 1/3 



(7) 



increases with r] while decreases with r]. Quite generally R^ < i?^, 
for T] > 10"* and vice versa for t] < 10"*. We should expect, therefore, that 
internal shocks will be important for heavily loaded "slow" fireballs while 
interaction with the ISM will be dominant for lightly loaded "fast" fireballs. 



3.3. BEAMING AND TIMING 

The fireball reaches relativistic velocities and the emitting source moves 
relativistically towards the observer. Thus, each observer detects radiation 
only from a narrow angle ~ I/7. This does not mean, however, that the 
overall GRB is beamed in a narrow angle. The angular spread of the GRB 
emission depends on the width of the emitting region which depending 
on the source model could be as large as iir (but not less than I/7). 

The duration of the burst depends on several factors. The fireball ap- 
pears as a narrow shell whose width could be as short at the original du- 
ration of the pulse at the source or significantly longer due to spreading of 
the pulse: 

^ 1 1.5 X 10-3sec Elilvi'^'n-^/^ for R > R^ 
[_ IQ-^^sec Rij otherwise 

A given observer will detect radiation from an angular scale 1 around his 
line of sight. This will lead to a typical duration of [36]: 

ATi ^ R-y/-/FC = 5sec ^- ji??;^^^^-^/^ 

The over all duration, AT = max(Ar;, AT^). The strong dependence of 
Ax on T] is an advantage, as it provides a possible explanation to the large 
variability in durations of GRBs. 



4. Binary Neutron Star Mergers 

Narayan, Piran and Shemi [37] and Phinney [38] estimated, using the 
three binary pulsars observed in the Galaxy, the rate of NS^Ms as ~ 
IQ-^'^^'^/year/galaxy. This estimate is probably too high, mostly due to 
the fact that the current estimates of the distance to Pulsars is larger than 
the one used in those calculation (see e.g. Bailes [39]). Still the rate is within 
the range of rates quoted earlier for GRBs (recall that the rates of long and 
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Figure 1. Logarithmic density contour lines at the end of the computation of the merger. 
The contours are logarithmic, at intervals of 0.25 dex 

short bursts are comparable). This is a crucial, albeit indirect, evidence for 
the NS^M model for GRBs. A small fraction of the total energy released 
in a NS^M suffices to power a GRB at cosmological distances. Thus the 
NS^Ms satisfies the two essential requirements from a viable model. 

There are still many open questions. The two basic ones are: (i) What is 
the energy transfer mechanism into the electromagnetic channel? (ii) Is the 
resulting fireball sufficiently free of baryons to reach the ultra-relativistic 
velocities (7>I0^) needed to produce a GRB? 

To address these issues we [5] developed a numerical code that follows 
neutron star binary mergers and calculates the thermodynamic conditions 
of the coalesced binary. The process of coalescence, from initial contact 
to the formation of an axially symmetric object, takes only a few orbital 
periods. Some material from the two neutron stars is shed, forming a thick 
disk around the central, coalesced object. The mass of this disk depends 
on the initial neutron star spins; higher spin rates resulting in greater mass 
loss, and thus more massive disks. For spin rates that are most likely to 
be applicable to real systems, the central coalesced object has a mass of 
2AMq, tantalizingly close to the maximum mass allowed by any neutron 
star equation of state for an object supported in part by rotation. Using 
a realistic nuclear equation of state we estimate the temperatures after 
the coalescence: the central object is at a temperature of ~ lOMeV, while 
shocks heat the disk to a temperature of 2-4MeV. 
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Fig. 1 depicts a typical density cut perpendicular to the equatorial plan. 
The disk is thick, almost toroidal; the material having expanded on heating 
through shocks. This disk surrounds a central object that is somewhat 
flattened due to its rapid rotation. An almost empty centrifugal funnel 
forms around the rotating axis and there is practically no material above the 
polar caps. This funnel provides a region in which a baryon free radiation- 
electron-position plasma could form [40]. Neutrinos and antineutrinos from 
the disk and from the polar caps would collide and annihilate preferentially 
in the funnel (the energy in the cm. frame is larger when the colliding v 
and V approach at obtuse angle, a condition that easily holds in the funnel). 
The numerical computations do not show any baryons in the funnels. The 
resolution of our computation is insufficient, however, to check whether the 
baryonic load in the funnel is low enough. The neutrinos radiation pressure 
on polar cap baryons can generate a baryonic wind that will load the flow. 
This effect depends strongly on the temperature of the polar caps [41, 42]. 
The estimated temperature from our computations, ~ 2MeV, is marginal. 
Our estimate is, however, least certain in low temperature regions like this. 

If the core does not collapse directly to a black hole, it will emit its 
thermal energy as neutrinos, vv e'^e~ could convert ~ 10~^ to 10""^ 
of the neutrino flux to electron-positron pairs and produce a GRB. The 
time scale for the neutrino burst is short enough to accommodate even the 
shortest rise times observed. Accretion of the disk onto the central object 
and magnetic held reconnection around the disk [31] are two additional 
energy sources that could power a GRB. This energy source can operate 
on a longer time scale independently of the dynamics of the central object. 

The numerical calculations support earlier suggestions [43] that the en- 
ergy release in anisotropic and that an empty funnel forms around the 
rotating axis of the binary system. The flreball is highly non spherical and 
it expands along the polar axis and forms a jet. This poses an immediate 
constraint on the model. If the width of the jet is 9 than we observe GRBs 
only from a fraction 20'^ of NS^Ms. The rates of GRBs and NS^Ms agree 
only if 9^0.2 unless the rate of NS^Ms is much higher than the current es- 
timates. In fact Tutukov and Yungelson [44] suggested recently that most 
neutron star binaries are born with a short orbital period and their life 
time is too short to detect them as binary pulsars. Such systems would 
escape the binary pulsar statistics and could increase by one or two orders 
of magnitude the rate of NS^M. 

The duration and spectra of GRBs vary from one burst to another. 
The flreball phase determines both the duration and the spectra of the 
bursts. The source might add variability by producing flreballs with differ- 
ent Lorentz factors and different initial durations. Within the funnel the 
baryonic load varies as a function of the angular position leading to varying 
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final Lorentz factors that, in turn, produce bursts with different durations 
and spectra. Another source of variability could arise from the interplay 
between the two energy sources in NS^Ms: Neutrino annihilation and ac- 
cretion energy of the disk. These mechanisms would operate on different 
time scale and produce different looking bursts. An additional source of 
diversity [5] is the distinction between systems that collapse directly to a 
black hole and those that undergo a longer rotating core phase. Finally, 
black hole-neutron star binaries are predicted to be as common as neutron 
star binaries[37]. Black hole neutron star mergers[45] would produce GRBs 
with different characteristics than NS^Ms. 

NS^M events can take place in a variety of host systems including dwarf 
galaxies, or even in the intergalactic space if in the process of formation 
the neutron star binary is ejected from the host galaxy [31]. Hence, unlike 
other cosmological models it is not essential that an optical counter parts 
will be observed in the location of GRBs [46]. A unique prediction of the 
NS^M model is that gravitational radiation signals from the final stages 
of the merger should accompany GRBs [18, 47] and vice versa (the latter 
is true only up to the anisotropic emission factor discussed earlier). This 
coincidence, which could serve to increase the sensitivity of the gravitational 
radiation detectors [47], would prove or disprove this model. 

I thank Ehud Cohen, Tsafrir KoUat, Ramesh Narayan and Eli Waxman 
for many helpful discussions. This research was supported by a BRF grant 
to the Hebrew University and by a NASA grant NAGS- 1904. 
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Paul C. Joss: I have a comment that is pertinent both to your talk and 
to that of the previous speaker. You have concentrated your attention on an 
essentially cosmological model for 7-ray bursts, and Gerry Fishman argued 
that there is no reason to invoke more than one mechanism for producing 
7-ray bursts and that, in particular. Galactic neutron stars are now ex- 
cluded as a source of 7-ray bursts. In fact, as shown by Gerry in his talk, 
7-ray bursts encompasses a very wide range of phenomenology in terms 
of overall bursts duration, temporal structure, spectral shape, and mean 
photon energy. I'd like to suggest that it would, if anything, be somewhat 
surprising if this broad range of phenomenology was due to a single un- 
derlying physical mechanism and that, in particular, it may be premature 
to exclude Galactic neutron star as the source of a modest fraction (per- 
haps 10%) of observed 7-ray bursts. The current situation is reminiscent of 
that surrounding cosmic x-ray bursts twenty years ago. THere with a much 
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narrower range of phenomenology, Occam's razor was frequently invoked 
to argue that all x-ray burst have a single physical source. This argument 
led theorists to speculate about a wide variety of exotic burst mechanisms, 
such as massive accreting black holes in the ores of Globular clusters, etc. 
When the work of Walter Lewin and his collaborators demonstrated that 
there were, in fact, two distinct types of x-ray bursts with different physical 
origins the theoretical situation was greatly clarified. Perhaps the same sort 
of clarification lies ahead in the future of our understanding of 7-ray bursts. 

Reply: I completely agree that the variability of GRBs suggests that 
there may be two distinct populations of sources. I tried to stress in my 
talk that the NS^M model allows for a lot of variability, but this might not 
be enough and one could easily hide a sub-class of about 10% of Galactic 
bursts amongst the observed GRBs. In fact I was led to the NS^M model 
while considering NS^Ms as gravitational radiation sources. I realized that 
NS^M release their binding energy in practically invisible channels and I 
was looking for a way to convert a fraction of this energy to electromagnetic 
energy. I was certain that if such a mechanism exists then NS^Ms would be 
observed via this channel. This has led me to the NS^M-GRB [7, 48] model 
in spite of the consensus that existed at the time that GRBs originate from 
Galactic neutron stars. In 1990 I was asked, in a seminar on this model at 
the lOA, if I really believed that GRBs are not Galactic. My reply was: 
"GRBs are such a diverse phenomenon that I could easily imagine that a 
small sub-class, say of 10%, are cosmological" . I am glad that now after 
the BATSE observations the situations has reversed and in reply to your 
question I would say that "GRBs are such a diverse phenomenon that I 
could easily imagine that a small sub-class, say of 10%, are Galactic". 

W. Kundt: You mentioned the difficulty for most models to account 
for the delayed super-high-energy detections (> lOGeV) from a few GRBs. 
In our contribution to the Untsville workshop, we explained them as the 
transient switch of a Geminga-like behaviour, cf. A&SS 200, 151 (1993). 

Reply: Several mechanisms, within the context of the fireball model, 
could produce the delayed super-high-energy detections (> lOGeV): (i) 
ultra-high energy photons produced at the shock of the fireball with the 
ISM (the delay is between the emission from the internal shocks and the 
ISM shock) [49]. (ii) Production of the ultra- high energy photons from 
interaction of extreme relativistic protons with a dense cloud of interstellar 
matter in the host galaxy (the delay arises from the fact that the trajectory 
of the protons is not directly towards us) [50]. (iii) Production of the ultra- 
high energy photons from interaction of extreme relativistic protons with 
intergalactic protons or CMBR photons. The time delay arises from a time 
of fiight delay of the protons. 
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